A Quantum Chemistry Plus Dynamical Mean-Field Approach for Correlated 

Insulators: Application to La2CuOi. 
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While the traditional local-density approximation (LDA) cannot describe Mott insulators, ab- 
initio determination of the Hubbard U , for example, limits LDA-plus dynamical mean field theory 
(DMFT) approaches. Here, we attempt to overcome these bottlenecks by achieving fusion of the 
quantum chemistry (QC) approach with DMFT. QC+DMFT supplants the LDA bandstructure by 
its QC counterpart as an input to DMFT. Using QC+DMFT, we show that undoped La2Cu04, is a 
d-Mott insulator, and qualitatively discuss the circulating current- and incoherent metal phase, at 
small but finite hole doping. Very good quantitative agreement with experimental photoemission- 
and optical spectra constitutes strong support for efficacy of QC-I-DMFT. Our work thus opens a 
new avenue for truly ab-initio correlation-based approaches to describe correlated electronic systems 
in general. 

PACS numbers: 71.28-|-d,71.30-|-h,72.10-d 
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High- Tc superconductivity (HTSC) in quasi-two di- 
mensional (2D) cuprates is an outstanding, unsolved 
problem in modern condensed matter physics. These 
materials are stochiometric Mott insulators (MI). Upon 
hole doping (x), a highly unusual, d-wave pseudogapped 
"metal" with strongly non-Fermi liquid (nFL) proper- 
ties smoothly evolves into the "strange metal" around 
optimal doping (xc), where singular responses reminis- 
cent of a IZ? Luttinger liquid are observed [l|. At lower 
T, c?-wave superconductivity, peaking around Xc, is seen. 
Upon overdoping, d-SC rapidly disappears along with a 
rapid crossover to a low-T FL metal. These unique ob- 
servations defy understanding in any FL picture, forcing 
one to search for non-FL alternatives. 

In reality, cuprates are charge-transfer (CT) "Mott" 
systems. (£) Strong particle-hole asymmetry is indeed 
necessary for a quantitative description. LDA based 
studies [3| suggest that the maximum T™°^ correlates 
with t' /t. Moreover, t" /t' is controlled by axial orbitals. 
(ii) The role of these apical orbitals in cuprates is ill- 
understood; they may be relevant for the "hidden order" 
in the pseudo-gap (PG) phase, presumed to be of the 
circulating current (CC) type 0, H, [1]. If true, how does 
this constrain a minimal model for cuprates 0, 0, S]? 
(Hi) A host of experiments [^, clearly reveal the 
k-space differentiation of quasiparticle (QP) states in 
doped cuprates. The normal-state PG has the same 
d-wave symmetry as the SC at lower T. Quantum 
oscillation measurements show that small hole pockets 
for small x evolve, possibly via multiple electron- and 
hole-like sheets 'll'|, into a full, Luttinger Fermi surface 
(FS) around optimal doping (xopt)- i'iv) Finally, near 
Xopt, singular low energy responses suggest a branch cut, 
rather than a pole-like analytic structure of the one- 
particle Green's function, G(k,Lu), near Ep [H- This 
goes hand- in-hand with FS reconstruction at Xopt 
Are these findings linked to a possible quantum critical 



point (QCP) at a; ~ Xopt, as in /-band compounds 13]? 
If so, does the PG state have a "hidden" order? De- 
veloping an ab-initio formulation capable of reconciling 
(iii) — (iv) above; i.e, the quantum oscillation (dHvA) 
data and angle-resolved photoemission (ARPES) disper- 
sion with one- and two particle dynamical responses is a 
challenge for theory and has hitherto been studied within 
effective models 14, HI]. If, however, the d-PG phase 
indeed carries CC order, [1, H, @] these issues must be 
studied using an extended Hubbard model involving both 
planar- and apical p — d states [1, Q . 

In light of (i) — (iv), inclusion of strong electronic 
correlations in a multi-band Hubbard model (with pla- 
nar and apical p — d states) is mandatory. Here, we 
extend earlier work [l6j to compute the dynamical re- 
sponses by marrying the quantum-chemical (QC) band 
structure with multi-orbital DMFT. We "derive" an ef- 
fective model for cuprates that, by construction, is con- 
sistent with (i) — (Hi) above. We show how both one- 
and two-particle spectra (PES and optics) in the CT- 
MI are quantitatively described by QC-I-DMFT. This is a 
novel theoretical route, hitherto unexplored, and avoids 
the "[/" problem in LDA-based approaches 17|, as de- 
tailed below. 

Earlier QC work is essentially a variational calcula- 
tion, using a multi-configurational wave-function (simi- 
lar to that used for the one-band Hubbard model [isj . 
but including local and nearest neighbor (n.n) p — d and 
d — d excitations), and quantitatively captures the strong 
renormalization of a doped carrier in a MI [l6l ]. Short 
range electronic (AF spin) correlations anisotropically 
renormalize the bare band structure, leading to k-space 
differentiation, as derived before in the one-band con- 
text [19I . [20| . and in good quantitative agreement with 
both, the ARPES dispersion [2l| , and the Shubnikov-de 
Haas data 11 , |23] ■ The nodal (N) QPs have predomi- 



nantly planar character, while the anti-nodal (AN) QPs 
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have significant mixing of apical — states. How- 
ever, the QC work cannot compute dynamical excitation 
spectra. To address (u) and {iv) above, the QC work 
needs to be "married" with DMFT/cluster-DMFT calcu- 
lations. Cluster-DMFT successfully reproduces various 
anomalous responses in cuprates, but within an effective 
model framework. Here, we harmonize the successes of 
the QC method and DMFT-like approaches for the CT- 
MI phase in an ab-initio framework; the doped case will 
be treated separately. 

Labelling the dj.2_y2 —pa and d^^ —pz bands found in 
the QC work by "1,2" leads to an effective "two-band" 
Hubbard model, H = Hq + Hi + Hmix, where Hq = 



i,a— 1,2 i.a^b i,a^b 

(1) 

and 



k,(T 



i„(k)(4i^Cfe2<T + h.c) 



(2) 



where tm(k) = tm{coskx-coaky) is the d-wave form factor 
associated with intersite, inter-orbital dx^-y^ — dz^ one- 
particle hopping. (Ai — A2) = 1.15 eV and U ~ 5.0 eV 
is chosen from the QC result {E^+i + Ejq^i ~ 2En — 
5.0 eV), while J ~ 0.6 eV is taken from atomic ta- 
bles, and U' ^ U - 2J = 3.8 eV. As for the hop- 
pings, t = 0.45 eV is chosen to be its "bare" value, while 
{f, t", = 0.01, 0.075, 0.2 eV) in H are renormahzed (by 
short-range static AF correlations) values taken from QC 
results [1^1 . This is because DMFT mainly renormalizes 
t, but the farther-neighbor hoppings are renormalized by 
non-local correlations beyond DMFT; we take the QC 
(static) renormalizations for these as an input into the 
DMFT machinery. Since QC already treats the effect 
of static correlations, [l3|, we do not include the static 
Hartree contribution in the DMFT treatment, thereby 
avoiding double counting of such terms. QC-I-DMFT 
now treats the effects of local dynamical correlations and 
non-local, static correlations on carrier dynamics in a sin- 
gle picture. The two dispersive bands are then given 
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'2t"{c2x + C2y) and 



with tz « t,t' ,t" ,tm- Here, 
x,y, z. While similar values for 



as ei(k) = —2t{cx + 
62 (k) = -2tz{cx + c 
Cq, =cos{ka) with a - 
the bare hoppings are found in LDA approaches, an ad- 
vantage of QC is an ab-initio estimate of the Hubbard U 
within a correlated formulation, avoiding the problems 
associated with LDA in this context [l3|- Also, in con- 
trast to LDA-based approaches, the QC work [16| gives, 
remarkably, a d-wave form-factor, and strong, local, d- 
shell correlations. Our two-band model is thus an ex- 
tended Anderson lattice model (EALM). Interestingly, 
Yin et al. Q derive a similar two-band-like model from 



LDA-\-U. In contrast to our QC results, however, with- 
out the d-wave hybridization. Within QC-I-DMFT, tm{k) 
gives a d-CT-Mott insulator, as we show below. 

We solve H within DMFT using the multi-orbital it- 
erated perturbation theory (MO-IPT) as the impurity 
solver [23|. Though not numerically "exact", it has 
many advantages: (i) it gives quantitatively accurate re- 
sults for band-fillings upto half-filling (23| at arbitrary 
T, (ii) self-energies can be easily extracted, and (iii) it 
is numerically very effcient. More "exact" solvers ei- 
ther cannot reach low T of interest (QMC) [13] or are 
prohibiti vely costly in real, multi-orbital cases (NRG,D- 
DMRG) [ii]. The relevant DMFT formahsm has been 
developed and used with very good success ^] for a va- 
riety of problems, and so we do not reproduce it here. 
The only input to the DMFT is the "free" density-of- 
states, given by p§(e) = J2k^(^ ~ ^a{kj)- We restrict 
ourselves to the quantum paramagnetic phase, above the 
3D Neel ordering temperature, T/v. With a non-local 
hybridization, the Green function (GF) is a 2 x 2 ma- 
trix, Gafc(k, Lj), (with a,b = 1,2) in orbital space, as are 
the /ocaZ self-energies, I]afc(k, w) = Sab(w). These equa- 
tions are similar to those appearing in the DMFT for 
the EALM The diagonal GFs (and self-energies) 

yield the total many body spectral function, A(k, uj) = 
(-l/7r)Im[l/(o.-S,(o.)-ea(k)-t2jk)Gb(k,a;))]. The lo- 
cal DOS is just Pq((^) = X)k ^a(}^, Using the spectral 
theorem, the off-diagonal spectral function, pi2(k,w) = 
(— l/7r)ImGi2(k, w), is easily seen to describe d-wave 
particle- hole (excitonic) order (cf. d-wave form of i,„(k)). 
The CT-MI (found below with QC-I-DMFT) thus has 
d-wave p-h order, shown by the fact that {c\^C2a) = 
(— I/tt) X)k / IniGi2(k, Ci;)dw = 0, as shown in the lower 
inset of Fig[l] . Our two-band Hamiltonian will also yield 
circulating current (CC) order at finite doping concen- 
tration, X, as proposed by Varma [27[. Here, however, 
the apical p — d link is crucial for CC order with finite 

A12 = i{icl„C2,^+e^_^,a " h.c)) = {Tf) ^ (sCC bclow) . 

This can be readily seen in the large-/7 limit of our model, 
where a second-order in t/U expansion gives the "ex- 
change" part as [13, [ill: 



H = 



U' 



<i,j> 



(1) 



(3) 



with P,^j = (S,.Sj- - 1/4), while tI^^ ^ [f^ 



aa ' "bb 
(2) 



KbKTrTJ - 1/4) + taahbiT+T- + h.c) and T^^^ 
tlb^T^T^ 



-TrTr) + {t,i,{taa-tbb)){TtT^+T^-T^^), where 



n2^)/2X 



CijC2, 



~ CjjCii. At mean- 
field level [27[ , (notice the difference in our T"s) this will 
yield a finite CC order, {Tf) ^ [13]. However, our 
CC pattern involves both planar and apical states, i.e., it 
is closer to that proposed by Weber et al. involving 
three oxygens on faces of the octahedra. While this sup- 
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ports the view 27| that p — d interactions are important 
in cuprates, our CC pattern is different in details. 
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FIG. 1: The "orbital resolved" and the total many-particle 
density of states (DOS) for the two-orbital model at half- 
filling, (n) — 3. Notice the appearance of "orbital" dependent 
Mott-Hubbard gaps, and the d-wave insulator, as explained 
in the text. 

We now present our results. In Fig[Tl we show the "un- 
perturbed" DOS (using QC) for our two-orbital model 
as dash-line curves. Clearly, the planar states dominate 
at Ep, which lies very close to a van-Hove singularity 
(vHs); the apical states lie much (~ 1.0 eV) lower. The 
many-body spectra within QC-I-DMFT, shown in bold 
lines, are dramatically different. Clear lower- and up- 
per Hubbard bands are visible in the DOS for both or- 
bitals. Dramatic and large-scale spectral weight transfer 
(SWT), driven by strong U, U', is readily manifest. The 
CT-Mott gap equals Amh = 1.1 eV. Our result can now 
be directly compared with the experimental photoemis- 
sion (PES) [2§| results for La2CuOi above Tjy. 

In Fig [51 we show this comparison. Quite remarkably, 
very good quantitative agreement with the PES spectrum 
upto —2.5 eV is clearly visible. The planar Zhang-Rice- 
like (ZR) states are the major contributor at low energy, 
but the apical states also contribute noticeably for en- 
ergies Lu > 1.3 eV. Further, from the self-energy (not 
shown), we estimate an effective mass enhancement of 
m* /mf, ~ 4.0, where uib is the LDA band mass. As a re- 
sult, we get {t,t',t") ~ (0.14,0.01,0.075) eV, and the re- 
sulting renormalized dispersion quantitatively fits the dis- 
persion of the one- hole ( "ZR-like" ) states in ARPES [16i] 
and the FS evolution with Ep shifts downward with 
hole doping. We show below that our QC-I-DMFT yields 
good quantitative agreement with the optical conductiv- 
ity as well. 

The theoretical optical conductivity is compared with 
the experimental data for La2Cu04 above T^. In our 
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FIG. 2: (Color Online) Comparison of the total one-hole spec- 
tral function obtained experimentally [2^ (triangles) with the 
QC-fDMFT spectrum with U = 5.0 eV (fuU line). Very good 
quantitative agreement is clearly seen up to binding energy 
~ 2.5 eV. 



two-band model, (t{oj) has two contributions: intraband 
transitions involving the planar- and apical states, and 
inter-hand contributions involving transitions betwen the 
two bands. The usual DMFT equation for now 
reflects both these processes J3(|, and reads 



depo{e) / dv p^{u; + iy)p^{v) 

(4) 

where p,{uj) = ELiPe„(w) = (-lA) ELi + 

fla-^a- Sa(w))]. 

Very good quantitative theory-experiment comparison 
is clearly seen in Fig [3l Specifically, the relatively sharp 
peak-like structure a.t H. — 2.0 eV, as well as the weaker 
shoulder (at fig ~ 1.2 eV) and the high-energy bump 
(at rjfc ~ 2.4 eV) are all in satisfying agreement with 
experiment Especially interesting is the shoulder- 

like feature at fig'- given that Amh = 1.1 eV, one would 
associate the 2.0 eV structure with the onset of the corre- 
sponding optical absorption feature. The shoulder at 
must then be interpreted in terms of a quasi-continuum 
"excitonic" feature pulled down below the Mott gap. In 
our two-band model, this arises directly from the inter- 
orbital transitions involving the planar- and apical bands. 
Thus, remarkably, QC+DMFT achieves very good quan- 
titative agreement with both, one- and two-particle dy- 
namical responses in the CT-MI phase, benchmarking 
its efficacy. 

What do we expect at finite doping, x7 For small x, 
we expect the lower- lying apical {d^^ — Pz) band to re- 
main Mott-localized, as is generic in MO-Hubbard mod- 
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els [20|, |23j, while the planar ZR band will selectively 
metallize. This would then be interpretable as an "or- 
bital selective" Mott transition (OSMT). One should 
then expect nodal fermionic QPs to dominate the re- 
sponses at small x [2^ , also found in earlier QC work [l^ . 
This OSMT would thus realize the famed N-AN di- 
chotomy Uc^ ubiquitous to cuprates. Eventually, 
around a critical doping, we expect the apical p — d band 
to metallize as well, as Ep shifts progressively downward 
with X. It is tempting to link the OSMT, where strong 
scattering between Mott localized (apical) and quasi- 
itinerant (planar) band carriers within DMFT would lead 
to low-energy infra-red singularities via the Anderson 
orthogonality catastrophe, to the doping driven avoided 
(pre-empted by d-SC) "quantum critical point" (QCP) 
around optimal doping, where low-energy singularities 
indeed dominate the "strange metal" More work 

is needed to check this theoretically, and to see whether 
this coincides with a T = melting of CC order 
These issues will be addressed in detail in forthcoming 
work. 
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FIG. 3: (Color Online) Comparison of the exp erimental 
optical conductivity (diamonds), taken from 31] with the 
QC-I-DMFT result (solid line) for La2CuOi. Very good quan- 
titative agreement, including description of the ~ 1.0 eV fea- 
ture below the charge-transfer-Mott gap, is clearly seen. 

In conclusion, we have proposed a new QC-fDMFT 
method to compute the correlated electronic structure, 
along with the one- and two-particle spectra, for Mott 
(CT-Mott) insulators. Extending our earlier QC re- 
sults, where very good agreement with the dispersion of 
one-hole "ZR-like" states and FS was found, we have 
shown how marrying QC with DMFT shows that, at 
"high" T > Tm, the CT-Mott insulator has d-wave or- 
der. The very good quantitative agreement with both 
PES and optical conductivity spectra in the insulating 



state of La2Cu04 constitutes strong support for this con- 
clusion. Our QC-I-DMFT modelling thus reconciles the 
ARPES and dHvA results with one- and two-particle 
spectral responses for the CT-MI phase. In light of 
recent ideas P, i, [13, [ill, these findings serve as an 
excellent starting point to study the physics of doped 
cuprates in detail within an ab-initio correlated approach. 
QC-I-DMFT should also serve as a new theoretical tool 
with wide application to other correlated systems of great 
current interest. 
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